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Abstract 

The transition form factors of the semileptonic decay Bg into scalar Dsq(2317) 
^ ' meson is calculated in the framework of 3-point QCD sum rule. The branching ratio 
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1 Introduction 

Semileptonic decays of mesons containing charm and beauty quarks constitute a very im- 
portant class of decays for determination of elements of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix, leptonic decay constants of heavy mesons, as well as for understanding the 
origin of CP violation which is related to the structure of the CKM matrix in the Standard 
Model (SM), because strong interactions involving semileptonic decays are more simple 
compared to that of hadronic decays. In semileptonic decays the long distance strong bind- 
ing dynamics can be parametrized as transition form factors, calculation of which is the 
main problem of these decays. For estimation of transition form factors some nonpertur- 
bative approach is needed. Several methods, such as the quark model, lattice QCD, QCD 
sum rules, large energy and effective heavy quark theories, have been used to calculate 
the transition form factors. Among these approaches the QCD sum rules occupy a special 
place, since it is based on the fundamental QCD Lagrangian. 

The QCD sum rules method [1] has been successfully applied to a wide variety of 
problems in hadron physics (see [2,3] and references therein). The semileptonic decay 
D -^ K^eUe is first studied in QCD sum rules with the three-point correlation function in 
[4]. Following this work, D+ -> K°e+z/e, D+ -^ K^*e+Ve [5], I> ^ Trez^e [6], £» -^ peve [7] 
and B -^ D{D*)('mi [8], D ^ (pli/i [9] are studied in the frame work of the same method. 
Note that Dg — * (j)iu£ decay is studied in light cone QCD sum rules [10], which is an 
alternative approach to the traditional QCD sum rules. 

Transition form factors appearing in semileptonic decays depend not only on the dy- 
namics of strong interactions between quarks in the initial and final state hadrons, but also 
on the structure of the hadrons involved in the semileptonic decays. Since in the present 
work we will consider a scalar meson in the final state, a few words about the scalar mesons 
are in order. The structure of the scalar mesons is still under debate. At present, there are 
different proposals about the nature of the scalar mesons that have been put forward, for ex- 
ample, their structures are considered as composed of qq, multiquark qqqq or meson-meson 
bound states. 

For studying the structure of scalar mesons much more experimental data and theoretical 
analysis are needed. The observation of two narrow resonances with charm and strangeness 
Ds,)(2317) in the Dgir invariant mass distribution [11-17] and Dsj{2460) in the D*7r'^ and 
Ds'~f mass distributions [12-14, 17-19] have raised discussion about the nature and quark 
content on these states [20]. The natural identification consists in considering these states 
as the scalar Ds(,(2317) and axial vector Dsj{2460) cs mesons, respectively. The result of 
the analysis of the radiative decays D,o(2317) -^ D^-f, D;^(2460) -^ D*j and D;^(2460) -^ 
Ds„ (2317)7 is in favor of the interpretation of quark content of these mesons as being 
ordinary cs mesons [21]. It is further observed that deviation of the transition amplitude 
in the infinite heavy quark limit compared to that of the finite quark mass case is quite 
sizeable. In the light of the Dg^ — » D*'-/ case, one can ask whether a change in the heavy 
quark effective theory results occurs for the Bg -^ D<j(,(2317)£z>£ decay for the finite quark 
mass, or not. In the present we try to answer this question. Note that this decay has been 
studied within the frame work of the heavy quark effective theory in [22]. 

In this work, we will study the semileptonic decay of Bg meson to D5q(2317) meson, i.e., 
Bg -^ Dsf^{2317)iue, in the frame work of QCD sum rules method. The paper is organized 



as follows: In section 2, we derive the sum rules for the transition form factors. Section 3 
is devoted to the numerical analysis and discussion, and contains a summary of our results 
and conclusions. 

2 Sum rules for the Bg -^ Dsq{2317) transition form 
factors 

The amplitude of the Bg -^ Dsi,{2317)£ui decay can be presented in the following form 

M = ^V,,va,{l - 75)^ {D,, |c7^(l - 75)61 Bs) , (1) 

where Vcb is the CKM matrix element which describes the transition of a 6 quark into a 
c quark, and for the sake of simplicity, in all further calculations we will denote Dsg(2317) as 
-Do- The main problem related to the calculation of the matrix element {Dq |c7^(1 — 75)6] B^). 
Obviously, the vector part of weak c^^{l — 75)6 current does not contribute to the matrix 
element considered above, which immediately follows from parity property of the hadrons 
and weak current; and only axial part of the weak current gives nonzero contribution. From 
the Lorentz invariance, this matrix element can be parametrized in terms of the form factors 
in the following way: 

(Do(p') |c7m75&| Bs{p)) = I [UV, + /_?J , (2) 

where ,f+{q^) and f-{q^) are the transition form factors, V^ = (p + p')m ^^^ % ~ (p^p')m- 
It is well known that /_ is proportional to the lepton mass, and especially for the r case it 
can substantially be important. For this reason we will take both of the form factors into 
consideration. In calculation of the form factors f+{q^) and f-{q^) we will employ the QCD 
sum rules method and proceed by considering the following correlator: 

Mp\p'\q')=i^ j d'xd^ye^^^'y-^^^^ {^\jDAy)Jt{^)U^)\^) , (3) 

where JDoiu) = sc, J5 = S756 and J^ = cj^'j^b are the interpolating currents of the scalar 
Do, Bs mesons and weak axial currents, respectively. 

Let us first calculate the phenomenological part of the correlator given in Eq. (3). 
This can be obtained by inserting the complete set of intermediate states with the same 
quantum number as the currents Jd^ and J5. isolating the pole terms of the lowest scalar 
and pseudoscalar Dq and Bg mesons, we get the following representation of the above- 
mentioned correlator 

n / 2 /2 2x (0 \Jdo\ Do) {Do 157^75^1 B^) {B^ \rs^,b\ 0) 

^ ^ (0 \Jd,\ h{p')) {h{p') |c7^756| Hjp)) {n{p) \s^,b\ 0) ^^^ 

The second term in Eq. (4) takes into account higher states and continuum contributions; 
and h and Ti form a complete set of mesons having the same quantum numbers as the 
ground state mesons. 



The matrix elements in Eq. (4) are defined in the standard way as: 

(0|JdoI^o) = foomDu , 

{B,\-s^,b\0) = -zJ^^^, (5) 

mb + rris 

where /d„ and fs^ are the leptonic decay constants of Dq and Bg mesons, respectively. 
Using (5), Eq. (4) can be written as 

^^(P^P'^ q'^) = - / ^' , ^% / 2 %T^ 2vf/+^M + /-^mI + excited states .(6) 

In accordance with the QCD sum rules philosophy, n^(p^,p'^, q'^) can be calculated from 
QCD side wit the help of the operator product expansion method (OPE) in the deep 
Euclidean region p"^ <^ {nii, + rric)'^ and p'"^ <^ {rric + rris)'^. Equating the two different 
representations of 11^ gives us sum rules for the form factors f+{q^) and f_{q^). 

The theoretical part of the correlator is calculated by means of OPE, and up to operators 
having dimension d=6, it is determined by the bare-loop and the power corrections (Figs. 
(l)-(3)), from the operators with d=3 {^pip)-, d=5 m^i^ipil)) and d=6 (i/jiIj)'^. Our calculation 
shows that d=4 operator (C^j,) gives very small contribution, and for this reason we do 
not consider it in the present work. 

In calculating the bare-loop contribution, we we first write the double dispersion rep- 
resentation as 

rper _ _ — ; ^g^gi ^ „' ' -— + subtraction terms . (7) 

(27r)2y i^s - p^){s' - p'^) ^ ' 

The spectral density Pi(s, s', g^) can be calculated from the usual Feynman integral with 
the help of Gutkovsky rule, i.e., by replacing the denominators of the propagators as follows: 

1 



p2 _ j^2 



-27ri6{p — m ) , 



which implies that all quarks are real. 

After standard calculations for the spectral densities we obtain: 

1 
+ 2ms{mb -mc){l + A + B)} 



p+{s, s', g2) = -— — — — (A' + A){1 + A + B) + [{m, + m,f - q^]{A + B) 



p„(s, s', q^) = ^^^^jj-—^{[A' + A + {mb + nicf - q^ + 2m, (m^ - m,)]iA - B) 

+ A'- A-2m,(mb + mc)} , (8) 



where A' = s — m^ and A = s — m^, and 






3 



and A(s, s', q^) = s^ + s'^ + q^ — 2sq^ — 2s' q^ — 2ss' . Here and in all following expressions, 
subscripts + and — correspond to the coefficients of the structures proportional to V^ and 
g^, respectively. Note that, in deriving Eqs. (7) and (8) we retain only linear terms in rus 
in order to take SU(3) violation effects into account, and higher order rus terms are all ne- 
glected. The integration region in for the perturbative contribution in Eq. (7) is determined 
from the condition that arguments of the three 5 functions might vanish simultaneously. 
The physical region in s and s' plane is described by the following inequalities: 

_ ^ 2ss' - 2mls +{s + s' - q^){ml - s) ^ 



According to the quark-hadron duality, the contribution of higher states in phenomenolog- 
ical part is parametrized in correspondence with the spectral density starting from s > sq 



and s' > s'q. 



1 mb-rric 



In what follows, we present the contributions of d=3, 5 and 6 operators. 
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rrib + rric 
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nicimh - mc)_ mh{mh - nic) 2 
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lyA fy-ii 






f^> = rnsirrib + rric) {ss) 



rrir 






12 






/+ = —^1^0 {ss) 



3ml{mi) — nic) 3ml{mb — rric) 2{mb — 2mc) 






/3 



(l^d /v*' 






/2 



2{2mb — rric) {fnb — iTLc){2ml + mbnic + 2ml ~ '^q'^) 



f- 



(5) 



12 



/y^ ry-i' 



ml {ss) 



ly* Zi lY*' ^ 



3m.l{m.b + m^c) 3m^(mf, + m.c) ^ 2{m.b + 3m.c) 






13 



/yO ly' 






/2 



2{m.c + 3m,b) {m,b + m.^) {2ml + '^h'n^c + 2m^ — 2q^) 



ry^ lyf 



ty^ fyi^ 



f- 



(6) 



81 



nag {ss)' 



12m?^{mb — m^c) 12m,l{m,b — m,c) 






/4 



ry^ryf 



4:m,c{m,b — m.c){2ml + m,bm,c + 2m.l — 2q^) 8m,c{7m.b — 8m.c 



/yZi /yi tj 



„/3 



Am^binib — m.c){2ml + m-bm^c + 2m^ — 2q^) 8m{,(8m6 — 7mc 



rytJry'^ 



ryO ry' 



A{bml - 20mbmc + 5ml + 2q^) 48 



ry ^ ly' ^ 
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48 

ryZiry' 



(9) 
(10) 

(11) 
(12) 



(13) 



(14) 



m^lm^s {ss) 



Qm^^ijUb — rTT-c) m^cirrib — m.c){5m.l + Am^bm.c + 5m^ — 5g^ 






/4 



ryZi ry'iJ 



2m,c{7m,b — ll«^c) rnbinib — mc){ml + Smbrric + m,l — q^) 






/3 



lyij ryi ^ 



5m^ — 20mf,mc + llm^ — 4g^ GmKnif, — rric) 2mh{5mh ~ ISnic 



(Y* ^ lY^' Zi 



ry*^ ly' 



lyO ly^f 



Ae) "i I- \2 



-12ml{mb + mc) 12ml{mb + rric) 






/4 



ly'i/yi' 



Amcirrib + mc)(2m^ + mbrric + 2m^ — 2g^) 8mc{7mb + 9mc) 



2t''3 









/3 



Ambimb + mc){2ml + mbrric + 2m1 — 2q^) 8mb{9mb + Tm^ 



/ytJ rf"' ^ 



rv-iry ff^' 



+ 



28(m^ - m^) 



/vi^ rt^i A 






np A rfTi 



GniKnib + rric) nicirnb + mc){5ml + Ambrric + 5ml — 5g^) 






/4 



/vi^ rff \) 



+ -minis {ss) 

y 

2m,c{7m.b + I2m.c) m^birrib + m.c)(m1 + Sm-bm^c + m.^ — q^) 






iz 



lyij qn' Zi 



m^ — ISrUbm^c — Trn^ 6m,l{m,b + nic) 4 



(>^^ lyf A 



lyiYTi A lyAfy-ii 



2m,b{12m,b + ISttt-^ 



(15) 



(16) 



where r = p^—m^, r' = p'^— m^. Note that, in the present work we neglect the 0{as) correc- 
tions to the bare loop. For consistency, we also neglect 0{as) corrections in determination 
of the leptonic decay constants fs^ and //)„ . 

Substitute Eqs. (8)-(16) into Eq. (7) and applying double Borel transformation B 
with respect to the variables p^ and p'^ {p^ -^ Mf, p'"^ -^ M2) in order to suppress the 
contributions of higher states and continuum, we get the following sum rules for the form 
factors /+ and /_: 



f±{g' 



{mb + mc 



{n.lJMl+mlJM^) 



-e^"'s. 



(2 



—^ / dsds'p±{s,s\q^)e 



(6) 



^(/f + /f + /f + /£ 



' „2x -s/M2-s'/M| 



(17) 



Double Borel transformation is implemented by the following expression: 



B- 



1 1 



\m+n 



1 1 

r(m) Tin) 



^-s/ml^-s' /ml 



3 Numerical analysis 

In this section we present our numerical analysis for the form factors f+{q^) and f-{q^). It 
follows from the expressions of these form factors that the main input parameters needed 
are the condensates, leptonic decay constants of Bg and Dq mesons, continuum thresholds 
So and s'q and Borel parameters Mf and M2. 

In further numerical analysis we choose the value of the condensates at a fixed renor- 
malization scale of about 1 GeV. The values of the condensates are taken from [22] and 



can be listed as follows: 



m 



fj,=l GeV 

{ss) 



-(240 ± 10 MeVf 

(0.8±0.1)(V'^) . 



The quark masses are taken to be rric = {n = rric) = 1.275 ± 0.015 GeV, ms{l GeV) ~ 
142 MeV [22], nib = (4.7 ± 0.1) GeV [23]. For the values of the leptonic decay constants 
of Bs and Dq mesons we use the results obtained from two-point QCD analysis: fs^ = 
209 ± 38 MeV [3] and foo = 225 ± 25 MeV [21]. The threshold parameters sq and s'q 
are also determined from the two point QCD sum rules: Sq = (35 ± 2) GeV^ [3] and 
s'q = (2.5 GeVy [21]. The Borel parameters Mf and M| are auxiliary quantities and 
therefore the results of physical quantities should not depend on them, if OPE can be 
calculated up to all order. In QCD sum rule method, OPE is truncated at some finite 
order. For this reason, we need to choose "working" regions for the Borel parameters where 
form factors are supposed to be practically independent of them. The choice of the working 
region for the Borel parameters Mf and M| should be based, on the one side, on the 
condition that that the continuum contribution should be small, and on the other side, the 
convergence of the power corrections. As a result of the above-mentioned conditions, the 
best stability is achieved for 10 GeV^ < Mf < 15 GeV^ and 4 GeV^ < M| < 7 GeV^. 

As a result of the above-summarized considerations, our analysis leads to th following 
predictions for the form factors aX q'^ = 0: 



u 
f- 



0.20 ±0.05 , 
-0.32 ±0.08 



:u 



The errors are can be attributed to the variation of the thresholds, decay constants, uncer- 
tainties in condensates and in quark masses. 

For completeness, we also present the results of HQET for the above-mentioned form 
factors that predicts, /+ = —0.37 and /_ = —0.15. Indeed, we observe that finite quark 
mass effects are essential in determination of the form factors. 

In order to estimate the width of B^ -^ Doiuf it is necessary to know the q"^ dependence 
of the form factors f+{q^) and f-{q^) in the whole physical region m^ < q^ < {ttib^ —'rrioay'- 
The q^ dependence of the form factors can be calculated from QCD sum rules (for details, see 
[5,6]). In the present work we have analyzed this dependence and used it in our numerical 
calculations. 

Having decided on the parametrization of the form factors, it is not difficult to obtain 
the expression for the differential decay rate 



dq^ 



1927r3 



m. 



-G'\V,b?\''\ 



m 



B.,' 



2 2\ 

mD„,(l ) 



nip 



q' 



X < - 



/-( 



(2g2 + mj 



+ 



U{q' 
{q^ ± 2m| 



(2m|^ ± 2m^^ - q 
Ml') 



2{ml-mljRe[U{q')r_{q')] 



2{mj,^-mj,Jq'Re[U{q')r_{q')] + f^q 




Taking into account the q^ dependence of the form factors /+ and /_, and performing 
integration over q^ and using the total hfe-time tb^ = 1-46 x 10~^^ s [24], we get for the 
branching ratio 

B{Bs -^ Doiue) ~ 10-=^ , (£ = e, ^) , 
B{Bs -^ Dqtu^) ~ 10"^ . 

Our predictions on the branching ratio is considerably smaller compared to that of the 
HQET results. 

At the end of our analysis, we would like to note that, although current B meson 
factories do not produce Bg meson, it is hoped to be possible to study the weak decays of 
Bg mesons in future-planned experiments at LHC The exclusive Bg — > DqIvi decay can 
be studied via the resonant production of the scalar Dq meson in the weak decay of the 
Bg meson. This analysis can give valuable essential information about the quark content 
of the scalar Dq meson. 

In summary, we study the semileptonic Bg -^ Dq£u£ decay in the framework of 3-point 
QCD sum rules. We calculate the transition form factors, and using these predictions, we 
estimate the branching ratio of the Bg -^ D^iui decay. 



References 

[1] M. A. Shifman, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B 147, 385 (1979). 

[2] L. J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rep. C 127, 1 (1985). 

[3] P. Colangelo, A. Khodjamirian, in "At the Frontier of Particle Physics/Handbook of 
QCD" ed. by M. Shifman (Worl Scientific, Singapore, 2001), Volume 3, 1495. 

[4] T. M. Ahev, V. L. Eletsky, Ya. I. Kogan, Sov. J. Nucl. Phys. 40, 527 (1884). 

[5] P. Ball, V. M. Braun, H. G. Dosch, Phys. Rev. D 44, 3567 (1991). 

[6] P. Ball, Phys. Rev. D 48, 3190 (1993). 

[7] A. A. Ovchinnikov and V. A. Slobodenyuk, Z. Phys. C 44, 433 (1989); V. N. Baiev, 
and A. Grozin, Z. Phys. C 47, 669 (1990). 

[8] A. A. Ovchinnikov, Sov. J. Nucl. Phys. 50, 519 (1889). 

[9] D. S. Du, J. W. Li, and M. Z. Yang, Eur.Phys. J. C 37, 173 (2004). 

[10] T. Ahev, A. Ozpineci, M. Savci, Eur.Phys. J. C 38, 85 (2004). 

[11] W. A. Bardeen, E. J. Eichten and C. T.HiU, Phys. Rev. D 68, 054024 (2003). 

[12] M. A. Nowak, M. Rho and I. Zahed, Acta Phys. Polon. B 35, 2377 (2004). 

[13] T. Barnes, F. E. Close and H. J. Lipkin, Phys. Rev. D 68, 054006 (2003). 

[14] E. van Beveren and G. Rupp, Phys. Rev. Lett. 91, 012003 (2003). 

[15] S. Godfrey, Phys. Lett. B 568, 254 (2003). 

[16] P. Colangelo and F. De Fazio, Phys. Lett. B 570, 180 (2003). 

[17] N. S. Craigie and J. Stern, Nucl. Phys. B 216, 209 (1983); V. M. Braun and I. E. 
Filyanov, Z. Phys. C 44, 157 (1989); V. L. Chernyak and L R. Zhitnitsky, Nucl. Phys. 
B 345 (1990) 137. 

[18] V. M. Belyaev, V. M. Braun, A. Khodjamirian and R. Ruckl, Phys. Rev. D 51, 6177 
(1995). 

[19] I. I. Balitsky, V. M. Braun and A. V. Kolesnichenko, Nucl. Phys. B 312, 509 (1989). 

[20] P. Colangelo, F. De Fazio and R. Ferrandes, Mod. Phys. Lett. A 19, 2083 (2004). 

[21] P. Colangelo, F. De Fazio, A. Ozpineci, Phys. Rev. D 72, 074004 (2005). 

[22] Ming Qui Huang, Phys. Rev. D 69, 114015 (2004). 

[23] B. L. loffe. Prog. Part. Nucl. Phys. 56, 232 (2006). 

[24] S. Eidelman et al, Particle Data Group, Phys. Lett. B 592,1 (2004). 



8 



